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ABSTRACT 


In the present work, finite element analysis of free 
vibrations of two rotationally periodic structures, namely, 
banded blades and finite continuous beam, has been considered, 
using the method of complex constraints. 

Mode shapes of rotationally periodic structures can be 
classified into the classes of symmetric modes, anti -symmetric 
modes and the rest. In the third category, the eigenvalues have 
multiplicity of two, i . e. , corresponding to every natural 
frequency, there exist an orthogonal pair of mode " shapes with 
eigenvector <u> and <u>. The complex eigenvector <Z> = <u>+i<u> 

is also an eigenvector of equation of motion. It is observed 

that, the deflection of one substructure has the same amplitude 
and a constant phase difference from the proceeding substructure. 
It is therefore possible to express the mass and stiffness 
matrices of complete structure, in terms of the first 
substructure only. Therefore a complete structure can be 

analysed by considering only the first substructure. This is 
called method of complex constraints. In this method, the 
elements of stiffness and mass matrices are complex numbers. A 
computer programme has been made to implement the method and it : 
has been illustrated with the help of two problems as mentioned 
above. IMSL subroutine DGVCCG has been used to find eigenvalues 
and associated eigenvectors of the complex generalized eigenvalue 
problem. 
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For the analysis of banded blades both axial and transverse 
vibrations have been considered, whereas for continuous beam only 
transverse vibrations have been considered. Results have been 
compared with the natural frequencies and mode shapes, obtained 
by standard finite element method. It is observed that, the 
present method gives quite accurate results, and computer memory 
requirements as well as computation time are much less as 
compared to conventional finite element method. 
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CHAPTER I 


INTRODUCTION 


1.1 GENERAL INTRODUCTION 

There are many instances in engineering wherein, the 
structures encountered are |>eriodic in nature. The most common 
example is a continuous beam with equal span length. Besides^ 
there are many structures which can be modelled as a continuous 
beam. Usually a simple model of a common aircraft fuselage is a 
beam on periodically spaced supports. Structures usually have 
two kinds of periodicities viz. linear and cyclic or rotational 
periodicity. In a periodic structure a single basic unit 
C substructure!) repeats itself within the structure. In a linear 
periodic structure the substructure repeats along a line. A 
cyclic periodic structure is a special case of periodic 
structures in which the starting point of the first substructure 
is the same as the last point of the last substructure. Linear 
and cyclic periodicity is discussed in more detail in subsequent 
chapters. Besides a continuous beam, another common example of a 
linear periodic structure is a plane truss. 

Examples of cyclic periodic structures occurring in nature 
are numerous. A simple model of banded blades of turbine discs 
and rotors of pumps and compressors, is also a cyclic periodic 
structure. Cooling towers with legs used in power plants is also 
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a cyclic periodic structure. All types of gears and rings can 
also be modelled as cyclic periodic structures, and in electrical 
industries the alternator end windings are also cyclic periodic 
in nature. 

The structures, mentioned so far, have periodicity properties 
in one direction only. There are some structures, which have 
periodicity properties in two directions. They are called doubly 
periodic structures. Grillage used for roofing, the ruber tube 
of a bicycle, are some of the existing doubly periodic 
structures. 

1,2 LITERATURE SURVEY 

The natural vibrations of various periodic structures have 
been studied by several authors. D. J. Mead £1] used wave 
propagation method to study the natural vibrations of infinite 
continuous beam. He described the phenomenon of wave propagation 
by a propagation constant. G. Sen Gupta £23 extended it to the 
case of finite continuous beam with various end conditions. They 
could reduce the free vibration analysis of complete structure to 
that of a substructure only, but their method is applicable only 
to the structures of simple geometry. B. L. Dhoopar , P. C. Gupta 
and B. P. Singh £33 used transfer matrix method to analyze 
vibrations of cable networks. They obtained the transfer matrix 
for individual members, from the closed form solution of the 
equations of motions. These transfer matrices were then used to 
derive the transfer matrix for the whole system. 
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The above menti oned methods have been empl oyed f or 
vibrations analysis of linear periodic structures. In 1979, 
D. L. Thomas £4] reduced the finite element analysis of 
vibrations of cyclic periodic structures, to that of one 
substructure only, by the. method of complex constraints. 
C. W. Cai , y. K. Cheung and H. C. Chan £63 proposed a 
transformation method for cyclic periodic structures. They 
expressed the response of a system in terms of a series of mode 
subspaces, multiplied by a transformation matrix, which they 
called the transformation matrix. For doubly cyclic periodic 

structures double ’U" transformation is required. The authors 
have employed transformation method to study the transverse 

vibrations of plane trussesC^^- 
1*3 SCOPE OF THE PRESENT WORK 

In the present work free vibration analysis of two cyclic 
periodic structures have been carried out. The first problem is 
the study of free vibrations of banded blades and the second 
problem deals with two span beam with various end conditions. 
Both the problems have been solved by two methods. Hie first 
method used is the method of complex constraints, proposed by D. 
L. Thomas. In the second method finite element assembly for 
complete structure is done to obtain the natural frequencies and 
associated mode shapes. 

1*3*1 Problem 1 : Banded blades 

The first problem analysed is a linear structure formed by 
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banded blades as shown in fig C3. laD. This structure can be used 
to model a specific stage of a turbine, with banded blades. If 
the radius of the turbine disc is made very large and the disc is 
considered to be rigid, then the structure can be idealized as a 
linear structure with fixed supports. 

The problem has been solved by the method of complex 
constraints as well as by the standard finite element method. 
Results are discussed in chapter III. 

Problem 2 t Finite continuous beam with three different end 
conditions 

The second problem deals with, natural frequencies and mode 
shapes of a two span beam with the following three sets of end 
condi ti ons: 

i3) Both extreme ends hinged, 
iiD Both extreme ends clamped, 
iiiD One end clamped, other hinged. 

For all of the above mentioned cases the beam does not 
satisfy the condition for cyclic periodicity. So the equivalent 
cyclic periodic structure is proposed, which takes care of all 
the three types of end conditions. It is shown in chapter III 
that for the above two-span beam the equivalent cyclic periodic 
structure is an eight span continuous beam resting on hinge 
supports. In this case also the problem has been solved by the 
method of complex constraints. And the results have been compared 
with the results obtained by finite element assembly for the 
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complete structure i , e. the two span beam. 

In both the cases the continuum has been discretized by beam 
elements. The finite element assembly for complete structure and 
subsequent application of boundary condition results in the 
following generalized eigenvalue equation, 

C tK] - CM3 D <u> = <0> Cl.l? 

where IKl and [M3 are stiffness and mass matrices respectively 
and <u> is the displacement vector. Solution of equation Cl. I!) 
gives eigenvalues co and associated eigenvectors <u>. 

The method of complex constraints reduces equation Cl,iy to 
that involving stiffness and mass matrices of first substructure 
only, with stiffness and mass matrices having its elements as 
complex numbers. The method directly gives orthogonal pairs of 
eigenvectors, for multiple eigenvalues. 

In chapter II, a brief discussion on periodic structure is 
given, followed by the finite element formulation for natural 
vibration of structures. Finally in this chapter the method of 
complex constraints has been detailed out. The discussion and 
the analysis of results is presented in chapter III. Chapter IV 
summarizes the conclusions and suggestions for future work. 
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CHAPTER II 


THEORY AND PROBLEM FORMULATION 

In engineering practice, there are many structures 
encountered, which are periodic in nature. Rings, banded blades 
of turbine discs , fuselage of an aircraft, rotors of pumps and 
compressors are some of the commonly occurring periodic 
structures. 

Z.l PERIODIC STRUCTURES 

By periodic, it is generally meant that, something appearing 
at regular interval of time or over some interval of length in 
space. In mathematics, a periodic function is defined as, a 
function fCx), such that, it is defined for all x and there 
exists a P > O, such that? 

fCx+PD = fCxD. 

P is called the period of the function and the function 
repeats after intervals of P. Similarly a periodic structure can 
be defined as a structure, in which, the geometry, defined over 
some length of the domain, repeats over the entire domain. A 
periodic structure can thus be divided in to a number of 
identical substructures. Once the substructure has been defined, 


6 



the whole structure can be generated by repetition of the 
substructure. 

2*1*1 Linear periodic structxire 

In a linear periodic structure, the substructure repeats 
itself along a line in space. Continuous beam and plane truss 
are linear periodic structures. 

2*1*2 Rotationally periodic structures 

A rotationally or cyclic periodic structure is one, in which 
the starting point of the first substructure is same as the last 
point of the last substructure. 

The projection of the centre line of the geometry of a 
cyclic periodic structure, in a plane perpendicular to it*s axis, 
forms either a circle or a polygon. If the geometry has been 
defined at some angle O, then it is identical at ^ where 

is 2n/N, and n & N are integers. N is the number of repeating 
units CsubstructureD , and n is any integer less than N. It 
follows that, once the geometry has been defined over a sector 
from & to C9-^y/y f the remainder of the structure can be obtained 
by repeated rotations through y/. Rings, banded blades of turbine 
discs, rotors of pumps and compressors, all are rotationally 
per i odi c st r uc t ur es . 
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Some of the linear periodic structures can also be modelled 
as rotational ly periodic structures, by imposing certain 

constraints at the two extreme ends, such that these two ends 
mathematically represent the same point. These constraints are 
equality of slopes, displacements and bending moments at the two 
ext r erne ends . 

Doubly periodic structures 

Structures having periodicity properties in two directions 
are called doubly periodic structures. A doubly periodic 
structure may possess linear or cyclic periodicity in two 
directions or it may have a combination of the two. 

Doubly cyclic periodic structures are not very common. 
However, an infinitely long cylindrical shell with periodic 
stiffening by ribs in axial and circular directions, is a doubly 
cyclic periodic structure. Similarly rubber tube of a bicycle, 
which has a torous shape, is also a doubly cyclic periodic 
structure. Grillage used for roofing, is another example of a 
doubly periodic structure. This has linear periodicity in two 
directions. As mentioned earlier such structures can be analysed 
as cyclic periodic structure by imposing constraints at the ends. 
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2*2 FEM FORMULATION OF SIMPLE STRUCTURES 


Structures considered in the present work consist of 
straight bars. Vibrations of straight bars are governed by a 
second order partial differential equation for longitudinal 
vibrations, and by a fourth order partial differential equation 
for transverse vibrations. Structures have infinite degrees of 
freedom and a closed form solution can be found only for 
structures having simple geometry. For finding solutions in the 
case of complex geometries approximate methods are used. Finite 
element method is the most powerful among the existing 
appr oxi mate methods . 

In the finite element analysis of free vibrations of a 
structure, the continuum which has infinite degrees of freedom is 
discretized. Discretization converts it into a finite degrees of 
freedom system. The element stiffness and mass matrices are then 
obtained. Finally, these element stiffness and mass matrices are 
assembled in the global ones. The partial differential equation 
thus reduces to a set of simultaneous algebraic equations, which 
in matrix form can be written as, 

C [K3 - CM] D <u> = <0> C2’ 1) 

where, CK] and CM3 are global stiffness and mass matrices 

respectively. The solution of the above equation gives natural 



frequencies o> and eigen vectors <u>. 

To obtain global stiffness and mass matrices, the element 
stiffness and mass matrices are first found. In the present work 
structures have been discretized by beam elements. Both, 
longitudinal and transverse vibrations have been considered. 
Here, the analysis has been carried out for one element, which is 
true for all the elements. 


Longitudinal vibrations 

The equation of motion, for longitudinal vibrations of bars 
is 

pA - EA a^u/at^ =o c2.e:> 

wh©r©, p is the mass density of the material of the beam, A is 

area of cross section and E is Young’s modulus of elasticity. 

Assuming linear variation for axial displacement u, 

CeD 


u = < Nj Ng > f 1 


C2. 3:> 


where u^ and Ug are node displacements and and Ng are shape 
functions given as. 


N . = 1 -x/h 

1 

Ng = x/h 


ca. 45 


where, h is the length of the element. 
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Substituting displacement function C2. 33 in equationCS. 23 
and applying Gal er kin’s method^ the following equation is 
obtai ned: 

J Ni [ ay aye cea au^^^ yayo - pA a^u'^^ yat^i dx = o 

Integration by parts and further simplification leads to. 


- - . .Cne3^ . .Cne3 

[m]<u> +Ck]<u> : 



C2. 53 


where, the nodal forces are given by. 


1 


'"-EA 


l=<=0 1 

1^2 J 


L EA 


Ix-h J 


and element stiffness and mass matrices are given by the 
foil owi ng equati on : 

[k] = J EA <N'> <N'>^dx 

Cm] = J pA <N> <N>'^ dx C2. 6D 

On substituting the expressions for shape functions and 
integrating equation C2. 6I>, the stiffness and mass matrices come 
out to be , 


Cm] = pAhyQ 


2 

1 


1 

2 


C2. 7aD 
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KA^h 


1 


-1 


Ek] = 


-1 1 


C2. 7bD 


2* 2« 2 Transverse vibrations of beams 

The equation of motion for transverse vibration of 
beams i s , 

El + pA d^v/dt^= O C2. 8D 

where, I is the area moment of inertia of cross section of beam, 
about an axis perpendicular to both axial and transverse 
directions. 

For a beam element, every node has two degrees of freedom 
viz. v and v" . 

The displacement function assumed is, 

Vj- *N3 Vg Vg' CE. 95 

where v^ , v^ ' and v^, v^' are transverse displacements and slopes 
at node 1 and node 2 respectively. 

Shape functions , and are given by 

P P 3 2 

N, =1-3 X /h + 2 X /h 
1 

Np = X- 2 x^/h 4- x^/h^ 

= 3 2 x^/h^ 

N = - x^/h^+ x^/h^ C2.10D 

4 

Substituting the assumed displacement function C2. 9!) in 
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equation of motion C2.8D and applying Gal er kin's method the 
f ol 1 owi ng equati on i s obtai ned , 

CEI + pA ydt^Z dx = o 

Integration by parts and further simplification results 


in. 


r 1 ^ n 1 ^ .CneD_ , . ,Cne> 

[m3<v > +lk] <v> - { f y 

where the nodal force vector is given by. 


C2. 113 




-M. 




where, and S^, are shear force and bending 

moment at node 1 and node 2 respectively. The element stiffness 

and mass matrices are given by. 

Cm] = JpA <N> <N>'^ dx, and 

Ck] = JEI <N"> <N">^ dx Ca.l25 

On substituting the expressions for shape functions C2.103 
and integrating equations C2. 123 the stiffness and mass matrices 
obtained are the following. 


[ml 


pAh/420 


■ 156 

22h 

54 

-13h ■ 


22h 

4h^ 

13h 

-3h^ 


54 

13h 

156 

-22h 

C2. 131) 

-13h 

-3h^ 

-22h 

4h^ 
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'la 

6h 

-12 

6h “ 


6h 

4h^ 

-6h 

2h2 


-12 

-6h 

12 

-6h 

C2. 1415 

6h 

2h^ 

-6h 

4h^ 



2# 2* 3 Assembly 

After the element stiffness and mass matrices have been set 
up, the next step is assembling these matrices in global ones 
Corresponding to a component of element stiffness and mass 
matrices, its location in global stiffness and mass matrices is 
determined using the connectivity. This is done for every 
element and the contributions of all elements are added. 

As an example, consider the case of bending vibrations with 
element stiffness and mass matrices given by equations C2.132> and 
C2.14!). After the assembly, equation C2. 115 will look like. 


CM] <v > + EK] <v> = 




C2. 155 


In the above equation, the shear forces and bending moments 
at the intermediate nodes gets cancelled because of equal and 
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opposite contributions from the neighbouring elements. Thus in 

the unknowns in equation C2. 153 are node displacements and shear 

force and bending moment terms , and S and M . 

^ linn 


Boundary Conditions 

The typical sets of boundary conditions are the 
f ollowing. 

13 One end is fixed and other free Cfig 2.1a 3 

In this case boundary conditions will be 


Vl=0 

V, ' =0 

1 


S =0 
n 

M =0 
n 


23 Both ends hinged Cfig 2.1b 3 

Boundary conditions in this case are 


Vl=0 

M, =0 
1 


V = O 

n 

M = O 

n 


33 One end fixed, other hinged Cfig 2.1c3 

In this case boundary conditions are 


vi=0 

V, ' =0 

1 


V = O 

n 

M = O 
n 


As an illustration of the application of boundary 
conditions, consider the first case. The first two rows and 
columns from equation C2.153 are removed to apply the conditions. 
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v^=0 and v' =0. Equation C2. 15) will reduces to, 

[M3 -CV) + [K3 <v> = <0> C2. 16) 

Now putting^ 
v= V q expC i wt ) 

V = Vq expCiwt) 

Equation Cl 6) reduces to, 

C [K3 - CM] ) <v> =<0> C2.17) 

which is same as Equation Cl. ID. This eigen value equation can 
be solved to give natural frequencies co and eigen vectors <v>. 

Equat i on Cl , 1 D has been ar r i ved at by appl yi ng boundar y 
conditions for case 1. For any other set of end conditions also, 
similar equation will be obtained, except that, the size of 
stiffness and mass matrices will be different, depending upon the 
number of degrees of freedom suppressed to take care of end 
condi ti ons . 

The size of stiffness and mass matrices in equation C2.15D 
increases as the number of substructures increase. Therefore, 
when the no of substructures are large, solving the eigenvalue 
equation C2.17D requires large computer memory and computational 
effort. 

The case of periodic structures, all the substructures are 
idl^ntiaal. This property of periodic structures can be 
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exploited to reduce the size of stiffness and mass matrices. 


D, L. Thomas L41 


has shown that the vibration analysis of 


periodic structures can be reduced to that of a single 
substructure only. Since the present work is based on the work 


of D. L. Thomas [43, his theory for rotationally periodic structure 


is being reproduced hero. 


2*3 FEM FORMULATION OF ROTATIONALLY PERIODIC STRUCTURES 

Consider a rotationally periodic structure consisting of N 
identical substructures. Let each substructure have J degrees of 
freedom. The total number of degrees of freedom will thus be NJ. 
A typical eigen vector <u> of the whole structure can be 
partitioned into N subvectors each having J components. The 


vector <u> can be written as, 
<u> = < <u > <u > . . 


^ cny. .T 

<u > > 


C2. isy 


where <u> 


is the eigen vector associated with the j 


substructure. Mode shapes for rotationally periodic structures 
can be classified in the following three categories; 


Ca> In the first category are 


those modes, in whioh eaoK 


subst,ructur® has bh© sa«i© mod© shap© as i-ts n©ighbo«rs i.©. , 


Cj5 cj-n:> 

u = u 


, f or all j 


For such mod© shap©s an ©i 


g©nv©ctor for th© complete 
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structure can be written as 


CID 


ci:> 


ci:> 


<u> = <u u u u > C2. 19a3) 

Cb> In this class, each substructure has the same mode shapes as 
its neighbours but is vibrating in anti-phase with them. 
Therefore the mode shape for the whole structure can be expressed 
as , 

<u> = < u -u u “U > C2.19b:> 

Cc> All other possible modes fall in the third category for 
which, 

Cj+10 ^ CjD ^ Cj+ID ^ . 

u^u and u =«i-u , for all j 

Class CcD modes are associated with multiple eigen values. 
For such cases, associated with each eigenvalue are a number of 
eigenvectors lying in a 2-D subspace of n -dimensional space. Any 
eigen vector lying in the 2-D subspace can be represented as a 
linear combination of two orthogonal vectors, contained in the 
same subspace. 


2* 3* 1 Properties of eigenvectors 

In this section class CcD eigenvectors will be 

considered in detail. The eigenvectors are normalized such that 
T 

<u> <u> = 1. The eigenvector, obtained by rotating <u> through 
one substructure, is also an eigen vector with the same eigen 
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value as that of <u>. 


Let this vector be <u'>, which can be 


expressed as. 


. . CND. . ci:>_ ^ C2:). 

<u'> = ■{ <u > <u > <u >. 


<u > } 


ca.ao:) 


Note that <u'> may or may not be orthogonal to <u>. Let the 
vector orthogonal to <u> be denoted by <u>, which has the same 
frequency as <u>. The vector <u> is also normalized such that 

^ j ^ 

<u> <u> = 1. The vector <u'> can be expressed as a linear 

combination of the two orthogonal eigenvectors <u> and <u> , 


<u'> = a<u> + b<u> 


C2. aiD 


where a and b are scalars. 


The vector <u'> has been obtained by rotating <u> through 
one substructure. If <u> is rotated in the same direction, 
through one substructure, the vector obtained will be orthogonal 
to <u'>. Let this vector be denoted by <u'>. This is 
associated with the same eigen value as that of <u'>. This can 
also be written as a linear combination of the two orthogonal 
vector <u> and <u> , in the following way, 

<u'> = -b<u> + a <u> C2.a2:> 


It is to be noted that when <u> , <u> and <u'> are 


normalized, a and b are respectively cos^ and sin^. Therefore 
<u'> automatically gets normalized. 

Substituting <u'> from equation C2. 202) in equation C2.215and 
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•then comparing the left and the right hand sides the following 
relations will be obtained: 


^ CND 
<u 

il 

A 

- ci:>. 

<u > 


b 

,-ClD 

<u 

, ci:) 
<u 

> = a 

, C2:)- 
<u > 

+ 

b 

,-C2D 

<u 

. C2:> 
<u 

> = a 

, C3D. 
<u > 

+ 

b 

,-C3D 

<u 


and so on.. J C2.23D 

The above relations in general can be written as follows: 

= a + b <u ,for all j Ca. 24aD 

Similarly one can write, 

<u = -b + a <u ,for all J C2. 24b:> 


. b <u 


<u = -b + a <u , f or all j 


2# 3# 2 Complex eigenvectors 

The eigenvectors <u> and <u> as well as <u' > and <u' > are 
real, normalized, orthogonal pairs of eigenvectors and are 
associated with same eigen value. The linear combination of 
these vectors can be written more conveniently in the form 
complex eigen vectors <Z> and <Z'> as, 

<2:> = <u> + i <u> 

<2'> = <u'> + i <u'> C2.25:> 


In terms of sub vectors one can write, 

+ i <u 






C2. 26:> 
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from t-he above equation, 

+ i <G 

Substituting from equation C2. 24a5 
expressions for <u*"'^ and <u 

b=sinv', equation C2. 26!) reduces to. 

From the above equation, 

<Z > = e ^ <Z > 

. 3I> . 2i y/ . _C1 . 

<Z > = e ^ <Z > 

iCN-lDv^ 

<Z > = e ^ <Z > 

In general, the above relations can be written as, 

<Z^'^^> = 2 < j < N+1 C2. 27bD 

if j=N+l, 

But, <Z*^'^^> = <z‘'^^>, therefore, 

iHw . iSnrr 

e ^ = 1 = e 

from which, 

or = Srrn/H C2. 28]> 


and C2. 24bD the 
and with a=cos^ and 

C2. 27a3 
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2* 3# 3 Valoos at * n* 

The V', defined in equation C2. 28!), may have M possible 
values. The n depends on a particular mode shape under 
consideration and is a measure of periodicity of the mode. If 
n^O and H/n is an integer then mode shapes repeats after H/n 

substructures. For example if n=2, H=8, i.e. , K/n = 4 means, 

C5D CID 
u = u 

cg:> ca:> 

u = u 

C7:> C33 

u = u 

C8D C4:> 
u = u 

•As the no of independent value of n is N, can have the 

f ol 1 owi ng i ndependent val ues . 

For even values of N : 


-aTTCN/2-i:)/N, . . -47T/N, -an/N, 

0. 

Bny'N , 47T/N , 

. . . arrCN/a- 


For odd values of N : 





-27T/CN-1D/2N, . . ~2fT/N, O, 2rT/N, 


2rTC N- 



The negative values of 

¥ 

correspond 

to anti 

clockwise 

rotation and positive values 

to 

clockwise 

rotation. 

All the 


anticlockwise eigenvectors are orthogonal pairs clockwise ones. 
Once the clockwise vectors <Z> are known, anticlockwis vectors 
can be generated simply by taking complex conjugate of <Z> and 
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■thon gonora.'ting oigon vectors for complo'to s'tructure using <Z>. 
Eigen value equa-tion is -theref ore not required to be solved for 
negative values of n. 

will therefore have the following values. 

For N even 

O , Stt/H , 4fr/H , .......... n 

For N odd 

O, ErrN, 4n/H, arr/CN-lD/aH. 

2* 3* 4 Method of Complex Constraints 

The generalized eigenvalue equation for free vibrations of 
structures is equation Cl. ID i . e. , 

C CK] - w^CM3 3 <u> = <0>. 

The complex eigen vector <Z> = <u> + i <u> is also a 

solution of the above equation i . e. , 

C [K3 - W^CM] D <Z> = <0> C2. 395 

In this section, the above equation will be reduced to that 
involving stiffness and mass matrices of the first substructure 
onl y. 

In general , a substructure may be connected to a number of 
other substructures. This number can be two in the simplest case 
or it may include all other substructures. Here the case in 
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which oach subs*truc*tur© is connected bo only bwo neighbouring 

subsbruct-ur es , has been considered. The common node bebween 

Cj-ID and substructure is considered to be a part of 

substructure. Similarly the common node between and 

substructure belongs to Cj+13^^ substructure. The stiffness and 

mass matrices for each substructure can be partitioned into a 

number of submatrices. The substructure contributes only to 

4-1.1 » *^1 > With Similar expressions for mass 

matrices. Similarly substructure contributes to ^ , 

j-1 j-1 

j^Cj ^ and The assembled stiffness matrix of 

j-lj J J"“l J J 

t h t h 

Cj-ID and j subs-tructur© will look like. 



bwo 


The diagonal submat rices will thus have contributions from 
substructures and off diagonal submatrices will have 
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contributions from one substructure only. The stiffness matrix 


for the complete structure can be written as, 

.,C1D 


11 


11 


12 


"IN 


"21 


o 


-1 

J J-1 


J-i 
"j-1 j-1 

JJ JJ 


J-1 
^j-1 J 


K 


J+1 J 


K'^ 


jJ-*-l 

J 

J+lJ-^1 


,CND 

"N1 


0 


CN-i:)_^^CND 
nn-j ^nN *^NN j 


C2. 3013 


Similarly mass ma^trix for -the complete structure can be 
formed. 

The vector <Z> in terms of subvectors can be written as. 


CID 

<Z> = { <z > 






C2. 31 :> 

From equations C2. 27^, <Z> in terms of first the 
substructure can be written as, 

ciD. iv^,_ci:), 2iv',_ci:>- iCN-i:)v',.7Ci :)- .t 

<2>=-{<z > e ^<Z > e ^<Z >...e <Z >> 

C2. 32:3 

Consider equation C2,292> with £K3 and similarly I M3 
given by equation C2. 3!)3> and <Z> by equation C 2. 3223. The first 
row of equation C2. 292), can be written as. 


11 11 


2.„C12) CN23 
-w [M^^ 


,Ci:3 

"12 


CN23 

^IN ^ 


M. 


,C123 


12 


-o — 


Mih 




,C15^ 

,ci:3 


e ' Z 
iCN-l23v'.7Cl23 


=0 

C2. 332) 





On si mpl i f"! cat-i on , oquat^ion C3.33D can roducas "to. 


-C [ C ^ ^ +{<r^ 'i ^ ^ ^ ^ N*“l D V' ir^ -j 

<[CKii ^^11 ^12 ^IM ^ 


-o.^[CM5“*Mjf = (2'“>=0 


C2. 341) 


,C1D, 


^ . , i C N— 1 w 

Consi der © ^ , 

^iCN— _ ^iNv' 

—iip 

= e ^ C2. 35:) 

To reduce ■the equation C2. 34> to that invol'ving subvector 
<.Z'"^"^y only, it is required to find the relation between the 
submatrices of the first and the last substructure. Since all 

the substructures are identical, 

cm _ ci:> 

1 ^22 

cm _ C15 

IN ^21 

^,cn:) _ ^,ci:> 

N1 12 

Similar relations apply to mass matrices also. Using the 
above relations and equation C2. 353), equation C2.34:) reduces to. 


2,^^ci:> ^jci:>^^ iv/^ciD^ iCN-i:)^/ . 

-« ^22 ^ ® ^12 ® ^21 
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C2. 36:> 

Equation C2. 36!) is an equation involving stiffness and mass 
matrix of first substructure only. This equation can be solved 
to give complex eigen vector Eigen vectors for the 

complete structure can be generated in the following way: 



Using the equation C2.27bD and C2. 28!) and the 

©1 Qenvoc'tors Tor all “the subs'Lruc'tur^s are ob'tained. Thus "the 

ei gfenvector for complo'te s'truclure •CZ>, is ob'tained. From 

equat-ion C2. 25D the real part of this vector is <u> and imaginary 

part is <u>. The two orthogonal vectors are thus obtained. 

It may be seen that, both stiffness and mass matrices, 

appearing in equation C2.36:> are hermitian, i . e. , real parts are 

symmetric and imaginary parts are skew symmetric. It is a 

property of hermitian matrices that the eigen values are real but 

the eigen vectors are complex. The solution of equation C2. 363 

2 

will therefore give real values of o> and complex eigen vector 

(MJ ORTHONORMALISATION 

The eigen vectors can be made CM3 orthonormal by finding the 
T 

quantity a=<u> CM]<u> and then dividing the the eigenvectors by 
square root of this number. The a can be obtained using the mass 
matrix of the first substructure only. For finding a, the 
relationship between the submatrices of substructures and that 
between the subvectors can be used. For the cases in which each 
substructure is connected to only two substructure, the 
expression for a is derived as follows: 

Let = X + i y 

From equations C2. 27bD and C2. 28D, the eigen vector for 
complete structure <u> is given by, 

T 

<u>=<X Xcosv'-Ysinv/ Xcos2v'-Ysi n2v' . . . XcosCN-i:)V'-ysinCN-lDv'> 
Since all substructures are identical i.e. , 



and 


12 

23 

- 


34 

N1 

^21 

= 

32 

- 

— W - _ 


43 

. . . 

IN 


the expression for a - <u> [M]<u>, after simplification comes out 
to be 


«=E^<XcosCj-i:)V/ -Y sinCj-i:>>'^CMjJ^+M22^:)<XcosCj-lDv^-YsinCj-i:>V/> 


N-1 

N 

<Xcosky-Ysinkv/>^ ^ 12 ^ 


-"E 

k=0 

E 

j=2 

<XcosC j -1 +kD v^^YsinC j — 1 

N-1 

N 

CXcosky— Ysinky> 


+E 

1- — 

E 

•i —O 

<XcosCk-»"l -j;>v^~YsinCk-«“l 


C2. 373 


Similar expression is obtained for <u> also. 


2» 3. 5 CoTi 5 >uter time and storage 

It is interesting to compare the computer time and storage 

required in a solution of the eigenvalue problem of a structure 

by the two methods. To find the solution by standard finite 

element method the equation Cl. 13 is solved whereas while using 

the method of complex constraints the equation C2. 363 is solved. 

The major part of the solution is finding the eigenvalues and 

eigen vectors, once the stiffness and mass matrices have been set 

up. Most of the subroutines take the CPU time proportional to 

the cube of the number of degrees of freedom. To solve equation 

3 3 

Cl. 13 the computational time will be proportional to N J . 
While using the method of complex constraints only J degrees of 
freedom will be used, but multiplications involved are complex 
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and ©ach complex multiplication is equivalent to four real 

multiplications. Therefore, the time required to solve equation 

3 

CS. 363>, for on© value of n, will be proportional to 4J . To 
analyse a problem completely the equation C2.36> is required to 
be solved for CN/2 + i:> times for N even Cwhich is the worse 
caseD so the computational time will be proportional to CN/2 + 15 
4J^ or C2N + 45J^. The ratio of CPU time required to solve the 
equation Cl. 15 and that required for solving equation C2. 365 will 


be. 


N^J^/C2N + 45 = N^/C2N+45. 

For N=4 this ratio will be 5.33, i.e. the computation 

required while using standard finite element is 5.33 times 
than that taken by the method of complex constraints. 


■time 

more 



CHAPTER III 


RESULTS AND DISCUSSIONS 


As discussed in the earlier chapters, two problems have been 
solved to test the efficacy of the method of complex constraints. 
The first problem analyzed is that of free vibrations of a 
simplified version of a banded blade structure, whereas the 
second problem is that of two span beam with three different end 
conditions. Results obtained by the method of complex 

constraintsC method and those obtained by the finite element 

assembly for complete structure C method E!) have been compared. 

3« 1 Banded Blades 

Idealizations made for the banded blade problem have been 
discussed in the first chapter under the heading **Scope of the 
present work**. The simplified structure consists of equally 
spaced four cantilever beams which are connected as a 
chain Cfig 3.1a:). The structure is divided into four 

substr uctures. The number of substructures could have been 
chosen as two or three also but two substructures are too less to 
understand the potential and behaviour of rotationally periodic 
structures. Three substructure have not been chosen because 3 
being an odd number, antisymmetric modes can not be obtained. 
Four substructures give all the three classes of modes, discussed 
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in cha.p'ter II. Though i't would have been b©*t*t©r bo includ© mor© 
“than four substrrucbur© bub as ib happanad bo b© bh© firsb probl©m 
ib was d©cid©d bo resbricb bo four subsbrucbures. The 

subsbrucbur© selecbed is shown in fig.C3. lb!). The sbrucbur© is 
cyclic periodic i.e. bh© deflecbion, bh© slop© and bh© curvabur© 
are same ab bh© firsb poinb of bh© firsb subsbrucbur© and bhe 
lasb poinb of bh© fourbh subsbrucbur©. 

3.1^1 IMPLEMENTATION DETAILS 

While solving by bh© subsbrucburing mebhod Cbhe mebhod 
of complex consbrainbsD , bh© finib© ©lemenb assembly for bhe 
firsb subsbrucbur© only is required, whereas while solving by bhe 
second mebhod C mebhod 2!) bh© finib© ©lemenb assembly for complebe 
sbrucbur© i.e. all four subsbrucbures is required. For bhe plane 
vibrabions of bhe sbrucbur© and hence bhe subsbrucbur©, each node 
will have bhree degrees of freedom namely, u,v and B. The 
elemenbal sbiffness and mass mabrices are of bh© following bype. 



rck] 

a 

0 



0 

Ck]^ 

b 



'[ m] 

a 

0 

-1 


1 



0 

Cm], 

b 







where [k] and [m3 , which are ©lemenb sbiffness and mass 

a a 

mabrices respecbively for axial vibrabions, are given by ©quabion 

C2. 73) and Ck3, and Cm3, , which are ©lemenb sbiffness and mass 
b b 
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mat. rices Tor bending vibrations are given by equations C2. 13!) and 


C2. 14D. The local displacement vector is of the form, 

,T 


<u> = < u. 


V, 

2 1 


"2 > 


and "the global displacement vector is of" the following type, 

T 


1 1 


<u> 

The el ement 




“2 ''2 ®2 


u V e y 

n n n 


stiffness and mass matrices have been 


transformed accordingly before assembly. The transformed local 
displacement vector is defined by, 

= CT] 

where [ T] is given by, 
for horizontal elements. 


■ 1 O 
0 0 


0 0 0 0 " 

0 10 0 



1 

O 


o o 

0 o 


0 0 0 0 
10 0 0 
0 0 10 
0 0 0 1 


for vertical elements. 


r 0 -1 



o 

0 

O 

O 

O 


0 O O O ■ 

0 0-10 

0 0 0 0 

1 o o o 

0 10 0 

0 0 0 1 


The transformed element stiffness and 
given by. 


mass matrices are 
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[Jc] = tT]^[k3^®^[T]^ , and 
[m3 = CT3^Cm3^®^[T3^. 

Tho ma.'t^rlal proper ■ties and geometrry d©*tails are as follows: 

E = 20 0 x'lO^ N/m^ 

p = 7000 Kg/m^ 

^ = 1.0X\0^ l.=lm 8tl =0. 3m 

A h V 

3»1*2 Coinparison of results obtained by the two methods 

The number of substructures in this case is 4, which is an 
even number. As discussed in chapter II, for even values of H, the 
independent values of n vary from O to H/2. So for H = 4, the 

independent values of n will be O, 1 and 2. Equation C2. 36:> is 

therefore solved for these three values of n. Table 1 gives 
first few natural frequencies of each case. 

According to equation C2.28!), = 2TTn/N. 

For n=0 v/=0 and 

for n=2 y/=n. 

The cases of n=0 and n==2 therefore represent symmetric 
C cl ass CaD 5 and antisymmetric C cl ass CbD modes respectively. 
The n==l case contains class CcD modes i.e. orthogonal pairs. 
Natural frequencies corresponding to n=l case are therefore 
repeated eigen values. Table 2 gives the comparison of natural 
frequencies obtained by the two methods. 

From Tables 1 and 2, the following points are observed. 
ly The frequencies obtained corresponding to the cases of n=0, 1 
and 2 include all the frequencies obtained by method 2 
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ii!) Frequencies ob*tai ned by met-hod 2, which ma*tch wibh "those of" 
n=l case are repeated eigenvalues. 

The mode shapes > obtained by sustructuring method, are [Ml 
orthogonal but not CM] orthonormal. CM] orthonormality property 
is important f"rom modal analysis point of view. The eigen 
vectors have been made CM] orthonormal by finding the quantity 

j 

a=<u> CM]<u>, using the expression C2. 373 and then dividing the 
the eigenvectors by squre root of this number. 

Fi gur esC3. 2— 3 . 1 03 show first three modes of each case of n. 
Figures C3.11—3.233 show the first thirteen modes obtained by 
method 2. Mode shapes corresponding to the cases of n=0 and n=2 
matches with the corresponding modes obtained by method 2. The 
two orthogonal modes of the case of n=l not match with 
corresponding modes obtained by method 2. But this is the case 
of repeated eigen values and it is well known that in such cases 
eigen vectors can not be represented uniquely. The eigen vectors 
in such cases lie in a 2-D subspace of n -dimensional space and 
that it is always possible to select a set of two orthogonal 
vectors which span the 2-D subspace. It follows that any vector 
lying in the 2-D subspace can be represented as a linear 
combination of the two orthogonal vectors. The substructuring 
method for n=l case directly gives two CM] orthogonal vectors 
<u> and <u>. The first natural frequency of n=l case matches 
with the third and the fourth natural frequency of method 2. It 
is observed that the following of the two eigen vectors 
corresponding to the first eigen value of the case n=l match with 
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■Lh© 3rd &nd ©igen v©ct.ors obt-ained by met-bod 3: 

O, 004-6010 — 0. 37576G1 -CuJ- = 3rd ©ig©n v®cbor and 

O. 3748304 •CuiJ’ + 0. 0403301 -Cuy — 4bli ©igen vector. 

Similarly, linear combinations can be made for other cases also. 

3.2. Finite continuous beam 

The second problem attempted is 2-span beam with the 
following three set of end conditions. 

Cal) Both extreme ends hinged CTyp© HHHD . 

CbD Both extreme ends fixed CTyp© FHF5 . 

CcD One end fixed, other hinged CTyp© FHHD . 

3.2.1 Construction of equivalent cyclic periodic structure 
Case 1 Both ends hinged Ctype-HHH) 

Figures C3. 24C^5 , C3, 24b> and C3. 24c:? the beam HHH and it's 
first two possible modes. Consider the first mode. It can be 
seen that the deflection, the slope, and the curvature are same 
at the two extreme ends. The conditions for cyclic periodicity 
is thus satisfied. It is observed that all the higher order odd 
modes have similar end conditions and hence they also satisfy the 
conditions for cyclic periodicity. Now, consider the second 
mode. It can be seen that the slope at the two extreme ends is 
not same. So the condition for cyclic periodicity is not 
satisfied. Similar is the situation with all higher order even 
modes. The method of complex constraint can not therefore be 
applied to this beam CHHH>. However, it is possible to construct 
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an equivalent cyclic periodic structure. The method of forming 
such a structure is discussed below. 

Consider again the second mode. In order construct an 
equivalent cyclic periodic structure, it is required to extend 
the original two span beam by its anti -symmetric image 
Cfig 3-273, so that the continuity of slope and bending moment 
Ccurvature3 is maintained. It will result in a four span beam 
resting on four hinged supports and which is a cyclic periodic 
structure. In this cyclic periodic structure thus obtained, the 
the deflection, the slope and the curvature are the same at the 
two extreme ends. Thus , to get the natural frequencies and mode 
shapes of the two span beam resting on three hinge supports, the 
four span beam should be analyzed while using the method of 
compl ex const r ai nts . 

Case 2 Both ends fixed Ctype FHF> 

The first two modes of this case are shown in fig C3. 25b3and 
C3.2Sc3. For the first mode, the slope and the deflection are 
same, both of these quantities being zero, at the two extreme 
ends. But, the curvature at the two ends are not same. 

Therefore the condition for cyclic periodicity is not satisfied. 
Similar is the case for all other odd modes. The second mode 
however satisfies the condition of cyclic periodicity. It is 
noticed that all the even modes satisfy the condition of cyclic 
periodicity. To analyse the two-span beam with this kind of end 
conditions, using the method of complex constraints, here also 
the equivalent cyclic periodic structure can be constructed. To 
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obtain the equivalent cyclic periodic structure it is required to 
extend the original beam by its anti symmetric image Ctig 3. SSD. 
This will result in a four span beam. It should be noted here 
that in order to be cyclic periodic, the equivalent structure has 
to have hinged support at all the four places. However, because 
of symmetry, certain mode shapes will have zero slopes at two or 
all the four supports. Such modes will belong to the case under 
consideration. The remaining modes will belong to the case 1 
i . e. , the HHH beam. 

Case 3 One end fixed, other hinged Ctype— FHH) 

In this case, one end is fixed and other is hinged. For 
every mode the slope at the fixed end will be zero and that at 
the extreme hinged end will be a non zero quantity. Therefore 
the condition for cyclic periodicity will never be satisfied for 
any mode. Figures C3.26bD show the first possible mode of this 
case. Here also an equivalent cyclic periodic structure can be 
constructed. For this purpose the original beam is first 
extended by its anti symmetric image. This will result in a four 
span beam with end conditions similar to odd modes of both ends 
fixed CFHHD beam. If this intermediate structure is again 
extended by its anti symmetric image, the structure obtained will 
be a cyclic periodic structure consisting of eight spans Cfig 
3. 2Q:>. Again, it should be noted that for the resulting 
structure to be rotational ly periodic, all the supports should be 
hinged. As stated earlier because of symmetry the slope at the 
two, four or all the supports will have zero values. Some of 
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those modos will belong to this typo of beam. 

is thus observed that for the first two cases the 
equivalent structure consists of four spans and for the third 
case it is an eight span beam. But it is not required to analyze 
the four span and the eight span beams separately. If the eight 
span beam is analyzed using the method of complex constraints, it 
will give natural frequencies and mode shapes of all the three 
set of boundary conditions of the two span beam. 

To take care of all the three types of end conditions of 
the two span beam, the continuous beam with eight spans, resting 
on hinge supports only, has been analyzed. 


3,2*2 IiiQ>lementation details 

In this case only transverse vibrations have been 
considered. The elemental stiffness and mass matrices are given 
by equations C2.13D & C2.14D. The local displacement vector is of 

the form, 

= < Vj ©I Vg ©3 >’■ 

and the global displacement vector is of the following type, 

T 


<v> ■= < Vi Vg 63 


V e > 
n n 


Therefore no transformation is required. 

The material properties and geometry details are as follows. 


£= 200X10'' 
p = 7800 


N/m 


Kg/m' 


^ = 1 .0 X 10"“^ 
A 


m 


Length of each span = 1 m. 
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Na'tura.l Troquoncies and mode shapes for eight.— span cyclic 
periodic cont-inuous beam have been obt.ained by subsbrucl.tjring 
mebhodC met/hod of complex conslraint-s!) . The subst.r uc*t ure is shown 
in f i g. C 3. 30D . To compare bhe results, natural frequencies and 
mode shapes for t.he t.wo span beam have been obtained by complete 
FEM assembly for the two span beam. 

3# 2* 3 Comparison of results obtained by the two methods 

In this case, the number of substructures are eight i.e. H=8, 
so the independent values of n will be, 
n=0, 1 , 2, 3 St 4. 

Table 3 gives first few eigenvalues corresponding to 
various cases of n. Table 4 enlists first few eigenvalues of the 
two span beams for the three different types of end conditions. 
Table 5 gives the comparison of natural frequencies of the two 
span beam with both extreme ends hinged Ctype HHHD case and the 
corresponding frequencies of equivalent cyclic periodic 
structure, i . e. the eight span beam. Similarly Tables 6 and 7 
give the comparison of natural frequencies of the two span beam 
for the other two sets of end conditions Ctype FHF St FHH2> , and 
the corresponding frequencies of equivalent cyclic periodic 
structure. 

The following points are be observed from Tables 3-7, 

1. As anticipated, the five cases Cn=0, 1, 2, 3, 4!) of eight 

span beam include natural frequencies of all the three types of 
the two span beam. 


39 



22. THo casos of n— -O and n— -4 includ© "tha na'tural i^r@<quoncios of 
the first two types of the two span beam only, i . e. HHH and FHF. 
This is because of the symmetric nature of the symmetric nature 
of end conditions, for these two types of the two span beam. In 
particular , these are odd frequencies of HHH type and even 
frequencies of FHF type for which the two-span beam itself is 
cyclic periodic. 

3. The natural frequencies of the case of n=2 also happen to 
match with the frequencies of HHH and FHF cases only. The case 
include those frequencies which are common to the type FHF and 
HHH. In particular, these are even frequencies of HHH type and 
odd frequencies of FHF type for which the two-span beam needs to 
be extended to make it cyclic periodic. 

4. The cases of n=l and n=3 include all the natural frequencies 
of the third case of the two-span beam, i . e. , FHH type. 

Mode shapes obtained in this case also are CM3 orthogonal 
but not CM] orthonormal. They have been made CM3 orthonormal 
using equation C2. 373) . 

Mode shapes obtained for eight span beams have been 
compared with corresponding modes of the two span beam by 
extending the modes of the two span beam over eight spans. 
Figures C3. 31-3. 40D show the first two modes corresponding to the 
five cases of n, including orthogonal modes. Figures C3. 41-3. 433 
show first four modes of each types of the two - span beam. 
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Mod® shapes of the th® cases of n*0 and n*4 

From equation C2.285, v^=27Tn/N, which means, 
for n=0, v^=0 and 

for n=4, 

Therefore the cases of n=0 and n=4 represent symmetric 
C cl ass Ca^D and anti symmetricC class CbDD modes respectively. 
These modes match with corresponding modes of HHH and FHF types 
of the two — span beam. As stated earlier the modes of the cases 
of n=0 and n=4 match with those modes of HHH and FHF types, for 
which the two span beam itself is cyclic periodic. Therefore, to 
compare these modes it is not required to extend these mode 
shapes of the two span beam, over eight spans. 

Mode shapes for ns=2 

Mode shapes of the case of n=2 correspond to class Cc!> modes 
i . e. , they are associated with repeating eigen values. As stated 
earlier, these modes match with those modes of symmetric end 
conditions for which the two span beam had to be extended in 
antisymmetric fashion to make it cyclic periodic, i.e, even modes 
of HHH type and odd modes of FHF type. The reason for this is as 
foil ows . 

For n=2, Therefore, the phase difference of various 

substructures with the first substructure will be as follows. 
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The mode shapes of the two span beams have been extended 
over eight spans Cfig 3. to compare them with the 

corresponding modes of the eight-span beam. It is seen that the 
two orthogonal modes obtained by substructuring method do not 
match with corresponding extended modes of the two span beams. 
It is observed that the linear combination of the two orthogonal 
vectors, obtained by the method of complex constraints, match 
wxth corresponding extended vector of the two span beam. The 
following linear combination of the two orthogonal vectors 
corresponds to the first natural frequency of the case of n=2, 
matches with 2nd vector of HHH case and the first vector of FHF 
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case. 


0.005870 <u> +0.3371595 <u> = 3nd eigenvector of FHF beam 

Mode shapes of the cases of n*l and n*3 

The cases of n==l and n=3 include mode shapes of the third 
type of the two span beam i.e. FHH type. These modes also belong 
to the class c eigen vectors i.e. the case of multiple eigen 
values. Here also to compare the results, mode shapes of two 
span beam Ctype FHHD have been extended over eight spans Cfigures 
3. 45^ and 3. 40]), In this case also, the modes of the two span beam 
and eight span beam do not match. The linear combinations have 
been made corresponding to first eigen values of the cases of n=l 
and n=3, which match with the corresponding extended modes of the 
two span beams. These linear combinations are as follows: 
for n=l 

0.9303335 <u> + 0.178898 <u> = 3rd eigen vector of FHH beam 
for n=3 

0.53095 <u> +0.47386 <u> = 1st eigen vector of FHH beam 

3* 2* 4 Identification procedure 

It is clear from the discussion in the previous section 
that the equivalent eight -span beam gives the modes of all the 
three types of the two-span beam. Now there is a need to develop 
a procedure to identify that which modes of various cases of 
correspond to which modes of the three types oi the two— span 
beam,' without solving the two— span beam by the standard finite 


n 
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element method. The procedure is as follows. 

Modes of the cases of n*0 and n“4 

Careful observation of mode shapes of these two cases 
of n reveals that, 

1 . All the even modes of n=0 and odd modes of n==4 have non zero 
slope at each support. This situation is possible only with the 
first type of the two-span beam, i.e. , type HHH. Therefore they 
correspond to HHH type beam. 

2. All the odd modes of n=0 and even modes of n=4 have zero 
slope at each support. This situation is possible only with FHF 
type beam. There fore these modes correspond to the modes of FHF 
type beam. 


Modes of the cases of n*l and n=3 

From equation 2.28, v^=2Trn/M, which means, 
for n=l , y/=ny4 and 
for n=3, v^=3rT/4. 

It is easy to see that any substructure is neither in 
phase or anti —phase with either of its neighbour. This situation 
will occur only in the third type of beam i . e. , FHH type beam. 
Therefore all the modes of n=l and n=3 correspond to the modes of 
FHH type beam. 

Modes of the case of n*2 


For n=2. 

V/=ny2, 

which means. 


. .ci:) 
<u> = 

- <U> 

= <u> = - 

. .C73 
<u> 

<u> = 

, ,C4:) 
- <u> 

= <U> = - 

<u> 


and 
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This "type of sit-uat-ion is possible only wit-h *the bwo 
span beam having symmetric end conditions. So the case of n=a 
should include the modes of FHF and HHH type of beam only. Also, 
n=2 is the case of multiple eigen values, it should therefore 
include those modes for which frequencies of FHF and HHH are 
same. 

In view of the above discussion, the procedure for 
identifying the modes of HHH, FHF and FHH types of beam is as 
f ol 1 ows . 

iD All the odd modes of n=4, even modes of n=0 and modes of n=2 
correspond to the modes of HHH type beam. by arranging the 
natural frequencies of these modes in increasing order, all the 
modes of HHH type beam can be identified. 

ii3> All the even modes of n=4, odd modes of n=0 and modes of n=2 
correspond to the modes of FHF type beam. Again by arranging the 
frequencies corresponding to these modes, all the the modes of 
FHF type beam can be identified. 

iiiD The modes of n=l and n=3 include modes of FHK type beam 
only. So all the modes of FHH beam can be obtained by arranging 
the frequencies of n=l and n=3 cases in increasing order. 


45 



% — 


V, rO 

s„=o 



fig 2.1a : Boundary conditions for a beam with 
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fig 2.1b : Boundary conditions for a beam with 
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hinged 

1 

2\ 

Vi =0 

VyisO 

V,'=0 

rhn- 0 


fig 2.1c : Boundary conditions for a beam tuith 
one end fixed, other hinged 
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fig 3.1b The substructure selected in probi 


Fig 3.1a The Banded blade structure 
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Fig. (3.24): Beam and F tg.O.PB): Beam FHF 
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fig (3.27 )equiva]ent cyclic periodic structure 


original structure 



fig (3.20) equivalent cyclic periodic strucutre 
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fig (3.29 ) equivalent cyclic periodic structure 
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fig 3.30 substructure selected in problem 2 
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Fig. ) ; Second mode o-f beam 


HHH extended over eight spans 
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extended over eight spans 
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Fig. (3.46 ): Second mode of beam FHH extended over eight spans 



TABLE 1 


Natural froquencios w CRad/s!) of banded blades for various values 
of n 


n=0 


n=l 


n=2 


168. 27474 
1158 . 06045 
2881 . 24023 
3006. 34082 
5430. 23388 


1003. 55780 
2246. 85424 
3085. 65722 
5519 . 82421 
6213 . 25300 


915. 92370 
1841. 84130 
3121. 18310 
5780 . 78515 
5870 . 05654 



TABLE 2 


Comparison of natural frequencies co CRad/sD of banded blade 
problem obtained by the two methods 


Natural frequencies 
obtained by standard 
Finite element method 


mode 

no. 


Natural frequencies of 

equi val ent 

rotational ly 

periodic 

structure 

n=0 

n=l 

n=2 

mod© 

(a> 

no. 

mode 

o> 

no. 

mode 

w 

no. 


168. 27488 
915. 94793 


168. 27478] 


915. 92370 


1003. 59039 
1003. 59039 
1158. 11508 
1842. 00316 
2247. 14539 
2247. 14395 
2881 . 79444 
3007. 15951 
3086. 54604 
3086. 54604 
3122. 12299 


1003. 55780 


1158. 06945 


2 1841. 841 3C 


2246, 85424 


3 2881 . 24023 

4 3006. 34082 


3 3085. 65722 


3 3121.1831' 











TABLE 3 


Natural frequencies w CRad/sD of the 8-span equivalent structure 


for different 

values of n. 




n=0 

n=l 

n=2 

n=3 

n=4 

1132. 9165 

1008. 7534 

780. 7320 

583. 0348 

499. 7668 

1999. 0738 

2169. 4196 

2530. 0820 

2918. 3645 

3122. 9494 

6122. 3676 

5834. 2714 

5278. 9067 

4751 . 3247 

4497. 9755 

7996. 6835 

8333. 3876 

9027. 6347 

9749. 6728 

10121. 1153 


Natural frequencies 

end conditions. 

TABLE 4 

CoD of “the 2-span 

CO i n Rad/s 

beam with three different 

Beam - HHH 

Beam - FHH 

Beam - FHF 

499. 7669 

583. 3035 

780. 7315 

780. 7315 

1008. 8753 

1132. 9165 

1999. 0738 

2169. 4197 

2530. 0823 

2530. 0823 

2918. 3645 

3122. 9495 

4497. 9755 

4751 . 3245 

5278. 9068 


99 



TABLE 5 


Comparison of natural frequencies ZoD of HHH beam with those 
various cases of n 



















TABLE 6 


Comparison of natural frequencies of FHF beam with those 

various cases of n 


Beam FHF 


Mode w 

No. C Rad/s :> 


Mode 

No. 


1 


3 


It 

o 

n — 2 

n = 4 

CO 

C Rad/s3 

Mode 

No. 

w 

C Rad/s3> 

Mode 

No. 

CRad/s3 


1 

780. 7320 



1132. 9165 

2 

2530. 0820 

2 

3122. 94£ 


3 

5278. 9067 

j 



6122. 3676 

4 

9027. 6347 




4 


10121. 115 




















TABLE 7 


Comparison of natural frequencies C«> of the Beam FHH with those 
of various cases of n 


Beam FHH j 

n = 1 

n = 3 

Mode 


Mode 

<A> 

Mode 

w 

No. 

C Rad/s2> 

No. 

CRad/sD 

No. 

CRad./s> 

1 

' '"1 

583. 0351 



■ 

i 

583. 0348 

2 

1008.7533 

1 

1008.7534 



3 

2169. 4197 

2 

2169. 4196 

i 


4 

2918. 3645 



2 

i 

2918.3645 

5 

4751 . 3245 



3 

4751 . 3247 

6 

5834. 2715 

3 

5834. 2714 



7 

8333. 3877 

4 

8333. 3876 



8 

9749. 6726 



■ 

9749. 6728 


1 OS 
















CHAPTER IV 


CONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK 


Th© method of complex constraints is applicable for 
rotationally periodic structures. This method can be applied to 
some engineering problems having linear periodicity properties, 
by appropriate modification of boundary conditions. In the 
present work this method has been applied to the following two 
linear periodic structures, by making them linear rotationally 
per i odi c : 

1. Banded blade structure, 

2. Finite continuous beam. 

Results have also been obtained by standard finite element 
method. 

Based on the discussion on results presented in chapter III, 
the following conclusions can be drawn. 

iD It is possible to apply the method of complex constraint to 
linear periodic structures by analysing their equivalent 
rotationally periodic structures. 

iiD Th© method directly gives orthogonal pairs of eigenvectors 
for multiple eigenvalues. Th© two orthogonal vectors can be used 

t-o span the 2— D vector subspace. 

iii3 By this method, it is easier to identify symmetric and 
antisymmetric modes, as they correspond to some particular values 


of n. 


ivl> It is possible to identify a particular eigenvector of the 
original structure from the modes of it*s equivalent rotational ly 
periodic structure by the suggested identification procedure. 
v3 Since the analysis is reduced to that of a single 
substructure, the computer memory requirements are much less as 
compared to the standard finite element method. conversely, for 
the same memory consumption the method of complex constraints 
will give more accurate results. The computation time also 
reduces considerably. The time saving is particularly 

significant when the number of substructures are large. Transfer 
matrix method also reduces the matrix size but as it uses an 
iterative procedure to find natural frequencies, there is not 
much saving of computation time. 

Suggestions for future work 

1. It would be interesting to extend this anal>'sis for finding 
the response of the periodic structures to external excitations, 
a. The method can be extended to the analysis of doubly periodic 

structures. 

3. In this method, in order to reduce the the analysis to that 
of a single substructure, use has been made of the relations 
between the submatrices of one substructure and the other 
substructures. For periodic structures these matrices are 

identical. If the relationship between the submatrices of 
substructures of almost periodic structures can be obtained then 
perhaps this method can be extended to the analysis of almost 
periodic structures also. 
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